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Abstract 

A new approach to the treatment of magnetic fluctuations in thin films of type 
I superconductors is introduced. Results for the dependence of free energy, specific 
heat and order parameter profile on the film thickness near the equilibrium phase 
transition point are obtained and discussed. 

Recent studies p!J EJ E] have shown that the Halperin-Lubensky-Ma effect (HLM) 0, 
namely, the fluctuation-driven change of the order of superconducting phase transition 
in a zero magnetic field is much stronger in quasi- two-dimensional (quasi-2D) type I 
superconducting films than in bulk (3D) samples. This fact gives new opportunities for an 
experimental verification of the effect in suitably chosen samples of thin superconducting 
films. Numerical values of thermodynamic quantities as latent heat, specific heat capacity 
and order parameter jump, have been calculated theoretically with the aim to support 
future experiments |3j- 

HLM effect occurs as a result of the interaction between the superconducting order pa- 
rameter ip and the vector potential A(x) = 8A{x) corresponding to the fluctuation part 
of the magnetic induction B = 5B; note, that the mean value B of the magnetic in- 
duction is assumed to be zero. When this interaction is neglected the superconducting 
phase transition in a zero external magnetic field is of second order. Taking into account 
the magnetic fluctuations A(x) and their gauge-invariant interaction with ip leads to an 
effective Ginzburg-Landau (GL) free energy that describes a first order phase transition 
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at the phase transition point T c0 = T C (B = 0). In type I superconductors the effect is 
stronger than in type II superconductors and can be investigated theoretically with the 
help of a self-consistent (mean-field-like) approximation, in which the scalar field tp can 
be considered uniform, i.e., independent of the spatial vector x. Here we shall use this 
approximation; a more detailed justification of this approach can be found in Refs. P3E]. 
The coupling between a scalar field, like ip in superconductors, and a vector gauge field, 
like A, is present also in models that describe other physical systems, for example, in 
scalar electrodynamics, certain liquid crystals and early universe phase transitions; see, 
e.g., [5 for a recent review. The problem is relevant also to quantum phase transitions in 
superconductors IE] and other systems jH] . 

In this letter we investigate the dimensional (2D-3D) crossover of phase transition prop- 
erties of type I superconducting films. Our treatment is based on preceding investigations 
of HLM effect |TJ where we have derived an effective energy of the superconduct- 
ing film from the Ginzburg-Landau model. We consider the effective free energy |3] and 
related thermodynamic quantities in the framework of a reliable and relatively simple 
variant of the crossover theory, used in the study of disordered films [10J. Because of 
considerable differences, this method can be applied to type I superconducting films only 
after an essential modification [H]. Here we outline a new general way of calculation of 
the thermodynamic functions dependence on the thickness of superconducting films. The 
present results, as summarized and discussed at the end of this paper, confirm and extend 
the recent investigations PjJ E] of the phase transition properties in thin films of type I 
superconductors. Our notations follow those in Ref. [T2|. 

We consider the following effective free energy density f(ip) = F(ip)/V of a thin super- 
conducting film of thickness L and volume V = (LiL 2 L ) PflEj: 

f^) = a\i,\ 2 + h -\^ + k B TJ[p^)] , (1) 



where 
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J(p)= I £-kS(k,p);, (2) 
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In Eqs. (JU)-©, p(ip) = PolV'l 2 ) where p = (8ire 2 /mc 2 ); a = a Q {T — T c0 ) and b > are 
the usual Landau parameters. They are related to the zero temperature coherence length 
£o = (^A^maoTco) 1 / 2 , the zero-temperature critical magnetic field H c0 = a T c0 (4:7r/b) 1 ^ 2 , 
and the initial (unrenormalized) critical temperature T c0 jT^j. The third term in the 
r.h.s. of Eq. (JTJ) describes the effect of the magnetic fluctuations [HE]- In Eqs. 
the integral J(p) and the sum S(k, p) over the wave vector q = (k, k ) are truncated 
by the upper cutoffs A and Ao. The finite cutoff A is introduced for the wave number 
k = \k\ of the wave vector k = (ki,k 2 ) and A stands for k . The free energy (JTJ)- 
(J3J) has been derived [HIS] under the condition k = [A(T)/£(T)] l/\/2, i.e. for type 
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I superconductors. Here, k is the Ginzburg-Landau number and the coherence length £ 
is given by £(T) = (^o/v^o)' where t = [(T — T c0 )/T c0 ] <C 1. The London penetration 
depth A(T) is A(T) = (Ao/y^o); Ao = corresponds to T = 0. 

In preceding investigations PQI21E1 of thin (quasi-2D) film, i.e., films obeying the condition 
ao <C Lq < A -1 , where ao is the lattice constant, only the (ko = 0)-term in the sum S 
has been taken into account. This approximation leads to a simple dependence of the 
thermodynamic quantities on the thickness Lq, but is not convenient when L varies in 
broad limits. Here we shall apply a more general treatment which gives the opportunity for 
the description of 2D-3D crossover from quasi-2D to 3D-superconductors L ~ Lp, (P — 
1,2). The choice of the cutoff A = (vr/£ ) is consistent with the general requirement 
£o < A(T) for the validity of GL free energy for type I superconductors (see, e.g., Ref. [T2] ) 
and specific requirements |3] for the validity of the Landau expansion of the free energy 
CQ)~© fo r extremely thin films (L <C £o)- 

As our study is based on the quasimacroscopic GL approach the second cutoff A should 
be again related to £o rather than to the lattice constant ao, i.e. A ~ (l/£o) 5 which means 
that phenomena at distances shorter than £o are excluded from our consideration. We 
shall assume that the lowest possible value of Ao is (tt/£o), as is for A, but we shall keep 
in mind that both Ao and A can be extended to infinity provided the main contribution 
to the integral J(p) and the sum S come from the long wavelength limit (g£ -C 1). 

In a close vicinity of the phase transition point T c0 from normal (ip = 0) to Meissner 
state > 0) the parameter p ~ ip 2 is small and the main contribution to the free 
energy / will be given by the terms in S with small wave vectors < A. This allows an 
approximate but reliable treatment of the 2D-3D crossover by expanding the summation 
over ko in (j^J to infinity - Ao ~ oo; for a justification, see Ref. [TTH lllj. A variant of the 
theory when A is kept finite (A = A = 7r/£o) can also be developed but the results are 
too complicated ^H]- Performing the summation and the integration in Eqs. (0)-© we 
obtain J(p) = (A 2 /27rL )/(p), where 



The integral (j3J) has a logarithmic divergence that corresponds to the infinite contribution 
of magnetic fluctuations to the free energy of normal phase (T c0 > 0, </? = 0). Such type 
of divergence is a common property of a lot of phase transition models. In the present 
case, as is in other systems, this divergence is irrelevant, because the divergent term 
does not depend on the order parameter ip and the free energy f(ip) is defined as the 
difference between the total free energies of the superconducting and normal phases: 



Introducing a dimensionless order parameter = (ip/?p ), where ip = (a T c0 /bY^ 2 is the 




fW = (fs-f N ). 
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value of if) at T = 0, we obtain the free energy in the form 

b, 



H 2 
8tt 



2t cp 2 + -\<pr + 2(1 + to)CI(fjup d 



(5) 



with /(/xv? 2 ) given by Eq. (J3J), /x = (l/7r/t) 2 , A = 7r/£ , and 

27r 2 A; B T c0 
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cQ 



From the equation of state {df /dip = 0) we find two possible phases: <^ 
superconducting phase (<p > 0) defined by the equation 



where 



+ , 2,(1 + to)CL £ 2 
^ i^A 2 ^ m ' 



K(z)= j\ y CoM\^^) 



(6) 

and the 
(7) 



io Vv + z 

The analysis of the stability condition (d 2 f/dip 2 > 0) shows that the normal phase is a 
minimum of f{<f) for t > 0, whereas the superconducting phase is a minimum of f(<p) 
provided 

.2i.;..^ (g) 
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1 > -(l+t )CL An 2 K(fi<p 2 oJ , 



where 
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o 2/ + ^ 



coth (| L kyjy + z) 



+ 



L A 



2sh 2 (| LqAv^T^) 



(10) 



The entropy jump As = (AS/V) = [—df(ipo)/dT] per unit volume at the equilibrium 
point of the phase transition T c ^ T c q is obtained in the form 



As(T c 



cQ 



1 + 
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where y2 c0 = <po(T c ) is the jump of the dimensionless order parameter at T c . The second 
term in As can be neglected. In fact, taking into account the equation f[(po(T c )] = for 
the equilibrium phase transition point T c we obtain that |C/(<^o)/Vol ^ s approximately 
equal to \t C Q + y? 2 /2|, where cp 2 and the dimensionless shift of the transition temperature 
t c o = to(T c ) are expected to be much smaller than unity. The latent heat Q = T c As(T c ) 
and the jump of the specific heat capacity at T c , AC = T c (dAS/dT) can be easily 
calculated with the help of Eq. (JTTJ). For this purpose we need the function <po(T), which 
cannot be obtained analytically from Eq. (|7jl. 

Even the relatively simple 2D-3D crossover scheme outlined in this paper cannot be used 
for an analytical calculation of the superconducting order parameter <f(T) and related 
thermodynamic quantities. But the simplicity of the scheme makes it very convenient for 
a numerical calculation of experimentally observable quantities in thin superconducting 
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films. Our numerical treatment of Eqs. (jSJ) and (J7J) shows a remarkably good confirmation 
of the results for 3D [2j and quasi-2D PI El superconductors. Here we want to mention that 
the reliability of present crossover scheme can be tested also analytically by considering 
the cases: LoA ^> 1, which corresponds to the "3D limit", and LoA <C 1, which gives 
the "quasi-2D limit". After taking the quasi-2D limit in Eqs. (0) and (J7|), we obtain the 
preceding results for (po(T) and the free energy f(ip) reported in Refs. PQEJ. The results 
in the 3D limit of Eqs. (0) and (JJJ) should be compared with the results for the free energy 
and the equation of state obtained when the continuum limit (J^ fc — ► J.) in the sum © 
is taken, which means to extend the integration to Ao = oo. In fact, we have done the 
calculations in the geometry of an "infinite cylinder" that corresponds to our finite-size 
treatment. Alternatively, the 2D-3D crossover can be described in the usual "spherical" 
geometry, when the cutoffs A and A are equal both in the sum (jSJ) and in the respective 
integral, calculated in the continuum limit (L A — > oo) along the small dimension L . The 
equations for the free energy and the equation of state in this variant of the theory are 
much more complicated and less convenient for a numerical evaluation of experimentally 
observable thermodynamic quantities. 
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